Interaction-driven spin precession in quantum-dot spin valves 
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We analyze spin-dependent transport through spin valves composed of an interacting quantum 
dot coupled to two ferromagnetic leads. The spin on the quantum dot and the linear conductance 
as a function of the relative angle 9 of the leads' magnetization directions is derived to lowest order 
in the dot-lead coupling strength. Due to the applied bias voltage spin accumulates on the quantum 
dot, which for finite charging energy experiences a torque, resulting in spin precession. The latter 
leads to a non-trivial, interaction-dependent, ^-dependence of the conductance. In particular, we 
find that the spin-valve effect is reduced for all 9 ^ tt. 

PACS numbers: 72.25.Mk,73.63.Kv,85.75.-d,73.23.Hk 
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Introduction. - The field of spin- or magnetoelectron- 
ics has attracted much interest, for both its beautiful 
fundamental physics and its potential applications. One 
famous spin-dependent transport phenomenon is the tun- 
nel magnetoresistance (TMR) in a spin-valve geometry, 
in which two ferromagnetic metals are separated by an 
insulating layer serving as a tunnel barrier |]J . The trans- 
mission through the barrier decreases as the relative an- 
gle 9 between the magnetizations of two ferromagnets is 
increased from to tt. Within a single-particle picture 
the ^-dependent part of the transmission can be shown 
to be proportional to cos 9. 

Transport based on tunneling has also been extensively 
studied in nanostructured devices such as semiconductor 
quantum dots (QDs) or metallic single-electron transis- 
tors. Recently, magnetotransport through those devices 
has attracted much interest. This includes normal or 
ferromagnetic metallic islands coupled to ferromagnetic 
leads jj, |5| as well as spin-dependent transport from fer- 
romagnets through QDs 0, 0, @ • Precession of a sin- 
gle magnetic atom spin in an external magnetic field has 
been detected, but only in the power spectrum of the 
tunneling current fol. Il0| - 

In this Letter, we study the effect of strong Coulomb 
interaction in a single-level QD (or a magnetic impurity 
(lUl attached to ferromagnetic leads on the average dot 
spin and the linear conductance in the weak dot-lead cou- 
pling limit (r -C knT, where T is the intrinsic line width 
of the dot levels) |l2j • We find an interaction-driven spin 
precession, even in the absence of an external magnetic 
field. This spin precession is predicted to be clearly visi- 
ble in the linear conductance as a reduction of the spin- 
valve effect and a nontrivial ^-dependence. For any non- 
parallel configuration, transport is reduced as compared 
to the parallel one (spin- valve effect). In the absence 
of Coulomb interaction, the ^-dependence follows simply 
cos 9 0,13- The presence of a finite charging energy, how- 
ever, leads to a reduction of the spin-valve effect. This 
can be understood by the interplay of spin accumulation 
caused by the bias voltage and an interaction-dependent 



spin torque due to an effective exchange interaction be- 
tween the spin in the dot and the leads, which in turn 
generates spin precession, detectable in the conductance. 



FM 



FM 




FIG. 1: A quantum-dot spin valve. A quantum dot (QD) 
is connected to two ferromagnetic leads (FM). The coordi- 
nate system we use is shown on the right. The magnetization 
directions (arrows) enclose an angle 9. 

The model. — We consider a small QD with one spin- 
degenerate energy level e participating in transport. The 
left and right lead are magnetized along Al and (see 
Fig- HJ, with relative angle 9. The total Hamiltonian 
is H = H dot + H L + H n + Ht.l + #t,r- The first part, 
Hdot = e^o- c o- c c + C/nj?ij, describes the dot energy level 
plus the charging energy U for double occupation of the 
dot. The leads are modeled by H r = e kal k a r k with 
r = L,R. For simplicity, we assume them to be half 
metallic, i.e., only majority spins have a finite density 
of states. Tunneling between lead and dot is described 

by -f/r.L = t^2k ( a Lfe c + + h.cA, where c + is the Fermi 
operator for an electron on the QD with spin along fiL. 
It is convenient to quantize the dot spin along the z- 
direction in the coordinate system defined by e x = (Al + 
n R )/|n L + n R |, e y = (n L - n R )/|n L - n R |, and e 2 = 
(n R x fiL)/|n R x £ll|. The tunnel Hamiltonian, then, is 



Ht.l — 



e^a{ kC] 



h.c 



(1) 



and ifp,R is the same but with L — * R and 9 — > —9. 
Due to tunneling the dot level acquires a finite width 
T = 2Tr\t\ 2 N, where N = N L = N R is the density of 
states of the majority spins of the left and right lead. 

With this choice of the quantization axis, the model 
studied here appears similar to those for Aharonov-Bohm 
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interferometers which contain a single-level QD in each 
arm 0,0,0]- In the present case, the two "arms" are 
labeled by the spin a =f , J, along the z-direction, and the 
Aharonov-Bohm phase due to an enclosed magnetic flux 
corresponds to the angle 9 between the leads' magnetiza- 
tions. The limit U — ► oo of the present model is, in fact, 
equivalent to the one studied in Sec. IV. C. of Ref. H3, 

Linear conductance and average spin. - We make use 
of the analogy between the quantum-dot spin valve and 
the Aharonov-Bohm interferometer and express the cur- 
rent in terms of Green's functions of the dot electrons, as 
shown in Eq. (4.3) of Ref. Q Here, we arc interested in 
first-order transport in T, for which the linear-response 
conductance G hn = {dl / dV)\ simplifies to 0] 
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Here, /(w) is the Fermi function, G TO '(w) arc the 
Fourier transforms of the Green's functions G > , (t) — 
-i(c a (t)cl(0)), G< a ,(t) = i(4,(0)c ff (t)), and Gp is the 
usual retarded Green's function. Contributions involving 
Imffjl^w), G^(w), and G^(w) are accounted for in the 
prefactor 2. Since T already appears explicitly in front of 
the integral, all Green's functions are to be taken to ze- 
roth order in r. In this limit, we find — (l/7r)Im G 1 j c j t (aj) = 



'11 

U), G> M 



(P° + Pl)6(w - e) + (Pf + P*)8{u> ■ 
2mP^8(u -e-U), and G<(w) = 2mP}S(uj - e), where 
P*, = (|x')(x|) are elements of the stationary density 
matrix (to zeroth order in T) of the quantum dot sub- 
system, with x, x' — ( em pty dot), |, J, (singly-occupied 
dot), and d (doubly-occupied dot). 

The main task is now to determine these density ma- 
trix elements to zeroth order in T. They contain as well 
the information about the average occupation and spin 
on the QD. The diagonal matrix elements, P*, are noth- 
ing but the probabilities to find the QD in state x, Le., 
the dot is empty with probability Pg, singly occupied 
with Pi = P| + Pp and doubly occupied with Pf. A 
finite spin can only emerge for single occupancy. The av- 
erage spin S = (S x , S y ,S z ) is related to the off-diagonal 
matrix element Pj via S. 



s z = (i/2)(p; -p l> . 

It is remarkable that on the r.h.s of Eq. J2J deriva- 
tives of Green's function with respect to bias voltage V 
appear. As a consequence, not only the equilibrium den- 
sity matrix elements enter the linear conductance, but 
also linear corrections in V are involved. In equilibrium, 
V = 0, the density matrix is diagonal with Pg = 1/Z, 
P T T = P/ = exp(-/3e)/Z, Pf = exp[-/3(2e + U)]/Z, with 



R-6 JP^ ■ Sy 



ImPj, and 



Z= l + 2exp(-/?e)+exp[-/?(2e+t/)] and = l/(k B T). 
As a consequence, the average spin on the QD vanishes 
at V = . With applied bias voltage, though, a finite 
spin can accumulate, which yields a finite (dS/dV)\ v _ . 

We now determine the needed density matrix elements 
by using the real-time transport theory developed in 
Ref. B The starting point is the generalized station- 
ary master equation in Liouvillc space, 



Xi,X2 



pXi 



E*, 1 

x 2 



XI 
X2 



0, 



(3) 



where xi an d Xi label the QD states, and e Xl and e X2 are 
the corresponding energies. The matrix elements Pj* 1 of 
the density matrix are connected to each other in Eq. © 

by terms E x }' Xl which can be viewed as generalized tran- 

J X 2 ,X2 ° 

sition rates in Liouville space. They are defined as irre- 
ducible self-energy parts of the propagation in Liouvillc 
space and are represented as diagram blocks on a Keldysh 
contour. For a detailed derivation of this diagrammatic 
language, the generalized master equation, and the rules 
on how to calculate a diagram we refer to Ref. [la- 
in the following, we write P = P + P + . . . and 
E = E + E + . . ., where P and E denote the equilib- 
rium limit, and P and E are the linear corrections in V. 
Using the symmetry of the model, we find the relations 



P* =P| andE x x x 



vx ,x 
^x,x" 



and P* = -P| and 



vx ,x 
^x.x" 



, where x is obtained from \ by the 
J.. For transitions from diagonal states 



vx ,x 
^x.x" ' 

transformation j 
in Liouville space to diagonal ones we find E x ' x , = 0. 
Finally, we drop all E terms which connect states in Li- 
ouville space that are not compatible, at least to lowest 
order in T. It turns out that it is sufficient to specify 
Eq. © for xi = X2 =1 as well as for xi =1, X2 =T- For 
the linear correction in V we get 



= 



T T ■ T 
_i_p0vQ4 



pi 



v ^ ■ J- 



yT,l 
T.T 



pi 



4,1 
4,T 



4,1 
4,1 



(4) 



(5) 



We evaluate the all necessary diagrams E*r* 2 explicitly. 
Eventually, we find the solution 



-i i eV ■ 
T Ak B T 



Pi cos 2 a{6) sin — 



(6) 



and P| = — PJ = iPi"tana(#). We used the definition 
A 9 



>1 



tana(0) = 



where A 



1 



/( C )+/( C +C0 C ° S 2' 



iRe 



2tt 



(7) 
and 



^(x) denotes the digamma function. This means that 
the spin accumulated in the QD is (for eV <C fc^T) 

eV 



^-P cos a (6) sin-, 



(8) 
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with a(9) being the angle enclosed by the quantum-dot 
spin and the y-axis, tana(0) = S z /S y . 

With the result for pi we are able to obtain the linear 
conductance. It can be written in the compact form 



G 



lin 



G 



lin, max 



1 



; a(9) sin 2 



(9) 



This equation together with the condition for a(9), 
Eq. 10, is the central result of this paper. The conduc- 
tance is maximal for parallel magnetization, 9 = 0. Its 
value is G lin > max = (ire 2 /h)(T/k B T)[l - f(e+ U)]f(e)[l - 
f(e) + f(e + U)}/[f(e) + l-f(e+U)}. 

We can straightforwardly generalize our theory to 
allow for arbitrary spin polarization p = (iV roa j — 
N m i n )/(N ma j + N min ) in the leads. In this case, we 
get G lin /G lin ' max = 1 - p 2 cos 2 a p (9) sin 2 (6/2) with 
tana p (0) = p tan a v =\(9), 

A reduction of the spin- valve effect has also been found 
for hybrid systems of ferromagnets with Luttinger liquids 
[l9| or normal metals j2jj. Both its physical origin and 
the ^-dependence of the conductance are different from 
our proposal, though. 

Results and discussion. - It is interesting to analyze 
how Coulomb interaction in the dot affects the spin accu- 
mulation in our model. In the absence of charging energy, 
U = 0, we find A = and a = 0, and the accumulated 
spin is along the y-direction, i.e., along — n^. A fi- 
nite charging energy, however, yields a rotation of the 
quantum-dot spin within the y — z-plane by an angle a, 
accompanied by a reduction of the total accumulated spin 
by cos a. The origin of the torque responsible for this 
spin rotation is an interaction-dependent effective spin 
exchange of the quantum-dot spin with the lead spins. 
In the subspace of single dot occupation an effective spin 
Hamiltonian can be derived from the full model by means 
of a Schricffer- Wolff transformation. If we now employ a 
mean-field picture and replace the lead spin operators by 
their average value, we end up with the simple effective 
Hamiltonian H cS = AT cos(9/2)S x for the subspace un- 
der consideration. In this effective model, a spin in the 
y — z-plane experiences a torque and starts to precess, as 
described by classical Bloch equations. Together with the 
rateri^ + ri^ d = r[l-/(e) + /(e+E7)] to diminish the 
z-component of the quantum-dot spin by changing the 
dot state from single occupation to an empty or doubly- 
occupied dot via tunneling, we get the Bloch equation 
dS z /dt = AT cos(9/2)S y - fT^o + Ti^ d )S z , from which 
we can extract the angle a(9) of the rotated spin in the 
stationary limit dS z /dt = 0, and we recover Eq. Q. 

We emphasize that both the rates for changing the 
number of dot electrons and the spin precession are of 
the same order in T, which is why the the angle a is T- 
independent. The two types of processes correspond to 
two different kinds of diagrams E x , 1:Xl . In spin-precession 
terms, all four labels Xi;X'i 1X2^X2 represent single oc- 
cupation, | , J.. They describe first-order virtual charge 



fluctuations during which the spin is rotated. In con- 
trast, tunneling rates which change the number of dot 
electrons are described by diagrams with \i = Xi = 

or Xi = Xi = d or X2 = X2 = 0, or X2 = X2 = d - lt is 
crucial for a consistent theory of first-order transport to 
include both types of diagrams in Eqs. Q and (J5J. 




FIG. 2: Upper panel: Linear conductance (normalized by 
F/ksT and plotted in units of e 2 /h) as a function of level 
position e for five different angles 6. Middle panel: Derivative 
of accumulated spin S with respect to bias voltage V normal- 
ized by ksT. Lower panel: angle a between the quantum-dot 
spin and the y-axis. In all panels we have chosen the charging 
energy J7/fcsT = 10 and half-metallic leads. 

The linear conductance as a function of the level en- 
ergy e is plotted in the upper panel of Fig. [3 for inter- 
action strength U/ksT = 10 and different values of the 
angle 9. For parallel magnetization, 9 = 0, there are two 
conductance peaks located near e = and e = —U, re- 
spectively. With increasing angle 9, transport is more 
and more suppressed due to the spin- valve effect. How- 
ever, this suppression is not uniform, as would be in the 
absence of charging energy. In contrast, the spin- valve 
effect is less pronounced in the valley between the two 
peaks, where the dot is dominantly singly occupied, and 
spin accumulation can occur. As a consequence, the two 
peaks move towards each other with increasing 9. 

The differential spin accumulation dS/dV in units of 
ksT is illustrated in the middle panel of Fig. [21 It is 
clear that single occupation of the dot is required for 
spin accumulation, i.e., the plotted signal is high in the 
valley between the two conductance peaks. 

As explained above, an effective exchange interaction 
between quantum-dot spin and spin of the leads yield a 
rotation of the accumulated spin in the y — z-plane by an 
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angle a. The lower panel of Fig. [5] depicts the evolution 
of a as a function of the level energy e. This angle is 
large in the valley between the conductance peaks, get- 
ting close to ±7r/2. A special point is e = —U/2, at 
which, due to particle-hole symmetry, the effective ex- 
change interaction vanishes. As a consequence, a shows 
a sharp transition from positive to negative values, ac- 
companied with a peak in the accumulated spin. 
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FIG. 3: Upper panel: Linear conductance as a function of 8 
for U/ksT — 10 and four different values of the level posi- 
tion. Lower panel: The same but normalized to the maximal 
conductance for parallel magnetization. 

In the upper panel of Fig. [2] we show the linear con- 
ductance as a function of 9 for four values of e. For 
e/ksT = 3 and 1, the dot is predominantly empty, and 
the ^-dependence of the conductance is almost harmonic. 
For e/ksT = —1 and —3, however, the spin-valve ef- 
fect is strongly reduced, and conductance is enhanced, 
except in the regime close to antiparallel magnetization, 
8 = 7T. The enhancement of conductance is related to the 
fact that the spin precession reduces the angle between 
the accumulated spin and the magnetization direction of 
the drain electrode. This is even better illustrated in 
the lower panel of Fig. [21 which shows the same curves 
but normalized to conductance at parallel magnetization 
(9 = 0). For e/ksT = —3, the conductance stays al- 
most flat over a broad range, and then establishes the 
spin- valve effect only in a small region around 9 = n. 

Finally, we comment that a finite spin-flip relaxation 
time T s f will reduce the spin- valve effect and limit its 
observability to T > r^ 1 (21|- The main prediction of 
our theory, the deviation from the cos#-law, will not be 
affected by r s f, as long as a ^-dependence is visible. 

To summarize, there are two pronounced features in 
the linear conductance which proves the spin precession 



proposed in this Letter: (i) the shift of two adjacent con- 
ductance peaks towards each other with increasing angle 
8 (Fig. |2J) , and (ii) strong deviation from the cos 9 law for 
the spin- valve effect (Fig. 121 . 
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